Entanglement dynamics in open two-qubit systems via diffusive quantum trajectories 
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We use quantum diffusive trajectories to prove that the time evolution of two-qubit entangle- 
ment under spontaneous emission can be fully characterized by optimal continuous monitoring. We 
analytically determine this optimal unraveling and derive a deterministic evolution equation for 
the system's concurrence. Furthermore, we propose an experiment to monitor the entanglement 
dynamics in bipartite two-level systems and to determine the disentanglement time from a single 
trajectory. 
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The fragility of entanglement with respect to decoher- 
ence persists as one of the major obstacles faced by most 
implementations of quantum information protocols. The 
unavoidable coupling of quantum systems to their sur- 
roundings is responsible for the deterioration of this most 
important quantum information resource, in ways that 
still escape quantification, and of which we have only a 
limited understanding. A complete account of this degra- 
dation process, even for the most elementary systems, is 
hampered, to a large extent, by the optimization problem 
inherent to the characterization of mixed state entangle- 
ment [TH3]- This is one of the major current challenges 
for theory and experiment. 

Recently, new light was shed on open systems entan- 
glement dynamics with the use of quantum trajectories 
[31 [S]. Based on numerical results, it was conjectured 
[S] that an optimal continuous measurement prescrip- 
tion generates the - in the above sense - optimal pure 
state decomposition of the system state, such that the 
mean entanglement yields the correct mixed state en- 
tanglement, at all times. Such optimal unravelings may 
help to overcome the mentioned experimental and the- 
oretical shortcomings, since they allow for a restriction 
of the optimization space, and define continuous mon- 
itoring prescriptions for the experimental measurement 
of entanglement. 

In the present paper we show that the use of diffusive 
quantum trajectories allows for an analytical proof of the 
existence of optimal unravelings in the sense of [5] , for two 
qubits coupled to a zero temperature environment. Our 
results, however, go beyond this point, complementing 
the program proposed in 0J [S] in three further direc- 
tions. First, we deduce a deterministic differential equa- 
tion for the evolution of entanglement in open systems. 
To our knowledge, the first of its kind. Second, we show 



that a single quantum trajectory unambiguously deter- 
mines the disentanglement time [5UlO). i.e. the time at 
which entanglement disappears completely. Finally, we 
propose an optical experimental realization of the opti- 
mal measurement strategy, for the actual monitoring of 
entanglement under incoherent dynamics. 

Consider a quantum system interacting with its en- 
vironment. Initially prepared in a pure state |^'o), it 
will evolve into a mixed state pt- The temporal change 
of the entanglement inscribed into the system state can 
be inferred by evaluation of some entanglement measure 
M{pt) at all times t, with [2] 

^iP)^, irliTp^M{^,), (1) 
{p.,l*i>}^ 

and A^(^') one of the known entanglement measures 
for pure states [3]. The infimum is understood over 
all possible decompositions of p into pure states, p ~ 
J2iPi\'^i){'^i\^ with positive < 1 and J2iPi = 1- This 
optimization turns into a formidable numerical prob- 
lem for higher dimensional or multipartite systems. The 
drawbacks, however, are not only computational. Even 
in those cases where ([T]) can be analytically evaluated 
[llj . e.g., for a system composed of two qubits, no clear 
physical interpretation of this optimal decomposition ex- 
ists such as to guide our intuition for the physics of mixed 
state entanglement. We believe that our present quan- 
tum trajectory approach indicates a possible strategy to 
overcome both these computational and conceptual ob- 
stacles. 

We analyse open quantum systems described by a 
Lindblad master equation [T^] 

N 

P = -^[H,p]+J2^kP. (2) 



2 



Here, the first term on the right hand side accounts for 
the unitary dynamics generated by the system Hamilto- 
nian H, while the superoperators Ck describe the effects 
of the environment on the system. Their action on an 
arbitrary state p is 

CkP = JkpJl - ^{JlJkP + pjUk) , 

and the operators Jk - which we will lump together as a 
vector J hereafter - are determined by the specific cou- 
pling between system and environment. 

Instead of attempting a direct solution of ^ for the 
mixed state p, we follow an alternative path in the spirit 
of [S] : A continuous monitoring of the environment yields 
information on the system, and induces a stochastic 
state evolution |13H15j conditioned on the measurement 
record. The unconditional state p is recovered upon aver- 
aging the conditioned state |^'c) over many independent 
realizations of the quantum trajectories at some given 
time t. Hence, the master equation is unraveled into 
quantum trajectories, with different unravelings arising 
from different ways to measure the environment |13| . It 
is this last feature which is at the core of our work. In 
what follows, we rigorously show that an optimal mea- 
surement scheme can be found that generates physically 
motivated decompositions of p into pure states, with min- 
imal entanglement average equal to the system entangle- 
ment as defined by Eq. ([T]). 

Our focus is on diffusive unravelings, yielding a contin- 
uous evolution for \^c)- The corresponding Ito stochastic 
equation for the state increment is p!4HT7] 



rf$^(i)(J-(J)c)|*c>, 



(3) 



where all expectation values are taken with respect to 
the conditional state I'I'c)- The first term describes the 
deterministic dynamics of the system, while the second, 
proportional to the vector d^ — (d^i, . . . , d^^v)^ com- 
posed of infinitesimal complex Wiener increments [18], 
encompasses the stochastic nature of the evolution. The 
stochastic process d^ has vanishing ensemble average, 
E[d^] — 0, with correlations 



d|d^ 



[ dt , d^d$J = udt , 



(4) 



where I is the identity matrix and u a complex symmetric 
matrix subject to the condition ||m||2 < 1 for the matrix 
two-norm 22 . 

Associated with the quantum diffusive process ([3| is 
the measurement record which conditions the evolution 
of l^'c), and which can be written as a vector of complex 
functions 



where each vector component represents a specific detec- 
tion event. The conditional change in Ivfc) is obtained af- 
ter using this experimental record to substitute the noise 
term d^ in ([3|. 

Notice that with the parametrization Q, a diffusive 
unraveling defined by ^ is completely specified once the 
correlation matrix u is given. Different unravelings, even 
though equivalent and leading to the same unconditional 
solution of the master equation ([2| , generate ensembles 
of trajectories with distinct statistical properties. In par- 
ticular, the optimal unraveling we are seeking must gen- 
erate, at each time step, an optimal pure state decom- 
position of pt which minimizes the average entanglement 
. Consequently, the necessary condition for an op- 
timal unraveling is that it has to minimize the average 
differential _E[(iA^] of the entanglement measure. This is 
achieved by optimization of the parameters which define 
u. 

Let us illustrate the ideas above by a specific exam- 
ple. Consider two qubits each coupled to a private bath, 
which we assume to induce decoherence by spontaneous 
emission (zero temperature reservoir). The correspond- 
ing Lindblad operators are J = {y/ja'^^ y/jl(^a^_^^y , 
where is the deexcitation operator of the ith qubit 
and 7 the spontaneous decay rate, which we assumed 
equal for both qubits. 

We will employ concurrence c |llj to quantify the sys- 
tem entanglement so far denoted by M. Analogous re- 
sults are obtained for other measures such as, e.g., entan- 
glement of formation [51lll|. The concurrence of a general 
pure state |*) = Voo|00) + ^oi|01) + iPw\10) + Vii|ll) 
is c(*) = \{^*\ay (g) ay\^')\ = 2|i/-ioV'oi - V'ooV'iil M- 
Under the dynamics generated by Eq. ([s]), the change of 
concurrence reads 



dcif — 



-Ul2 



dt 



(6) 



2c*Re((J^>*d^) 



where = c(^'), and = (**|(Ty ® CTj^I^') [TT]- 

The optimal unraveling is now found by minimization 
of the ensemble average E[dc<};]. This is accomplished 
by imposing that, for all trajectories and all times, the 
deterministic term in the first line of (|6| be minimal. It 
is readily verified that this is achieved by setting 



,(opt) _ 







_g»fopt 





(7) 



Y"^ dt = {J^u + J^)c dt + d$J , 



(5) 



with ^opt = arg(c^'(/'ii)- This choice of parameters im- 
plies that the optimal unravelling, determined by the 
phase of Q^si = c^V'iij is time independent: The equation 
of motion of Q^j, in its Ito form can be derived as 

de* = [2j\iPiifui2 + e*Re(/^.)] dt + e^Re{gld$,) , 

where and are functions of j^*). Note that this 
equation depends on the unraveling through ui2. Upon 
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substitution of u(°p*\ it reduces to 

dQ,^ = e,^Rc{f^ dt + gld^) , 

with /^f a new function of the state, and dQq, thus has the 
same phase as 9*. Hence, along trajectories generated 
with the optimal unraveling, the phase of 0^ remains 
unchanged - its value is fixed by the initial state. This 
feature has at least two significant consequences: 

(a) First, the optimal decompositions for the uncondi- 
tional state pt are singled out as ensembles of pure 
states which all exhibit the same phase of 8*, sig- 
naling the existence of a subtle relation between the 
temporal entanglement evolution and the charac- 
teristics of the optimal decompositions in Eq. Q. 

(b) Second, it implies that the optimal unraveling is 
time independent, with 

6lopt = arg(e*J = arg(4^V?i(0)) 

determined by the initial state |^'o) only. This re- 
markable property, observed numerically in [5] and 
analytically confirmed here, is of utmost impor- 
tance for the design of experimental setups to di- 
rectly monitor the entanglement evolution (see be- 
low): the monitoring strategy associated with the 
optimal unraveling is settled once the initial state 
is chosen, and does not require any steering along 
with the system's temporal evolution. 

As a last step, we still need to show that the optimal 
unraveling defined above indeed coincides with the state's 
mixed state entanglement as defined in ([T]), since the set 
of physically realizable decompositions of p is only a sub- 
set of all possible decompositions [H]. Substitution of 
u^i^^"^ into Q and a subsequent ensemble average leads 
to the equation of motion 



dE[c^] 
dt 



(8) 



The evolution equation of .©[l^/iup] is unrav- 
eling independent and can be derived to read 
dE[\^pii\'^]/dt = -27£:[|i/'iiP], which integrates to 
^[IV'iiP] = |V'ii(0)pe-27t. Thus Eq. ^ reduces to 

^^ = -7i?M-27|^n(0)|V2-*, (9) 

a deterministic equation of motion for the average con- 
currence, with parameters that depend on j'i'o) only. The 
solution of this equation is 

E[c^]{t, ^o) - e-^* [c(0, vPo) - 21^^11(0)12(1 - e"^*)] , 

(10) 

which coincides with the known result |SJ [50] for the spe- 
cific physical scenario under study here. 
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FIG. 1: Time evolution of concurrence for initial states j^'o) = 
(iVSlOO) - jOl) + j|10) 4- |ll))/\/8 (dashed line) and jl-o) = 
(|00) - jOl) -f ijlO) +iV5|ll))/v^ (solid line) under incoherent 
coupling to a zero temperature environment (time is rescaled 
as p = 1 — e"''*). Both states have the same initial concurrence 
c(^'o) = 0.809. The average concurrence evaluated using 500 
trajectories with an optimal unraveling faithfully reproduces 
the exact result [cf. Eq. (lOl]. Irregular lines represent the 
pure stateconcurrence along individual quantum trajectories 
for each initial state. The inset shows the time evolution of 
6 — arg(cji/)ii) along these trajectories. 



In addition to this coincidence between the average en- 
tanglement generated by optimal unraveling and mixed 
state entanglement according to ([l]), even single quantum 
state trajectories carry relevant information on the dis- 
entanglement process induced by the environment cou- 



pling: As evident from Eq. (10), entangled initial states 
with c(0, ^o) < 2|'(/'ii(0)P turn separable at finite times 

< oo, while they only asymptotically reach separabil- 
ity otherwise 121120] (see Fig.[T]). Since the concurrence 
of a two-qubit state is given by c{p) = max{0,A(p)} 
[llj . where A(p) — — X)i=2 a-^d \i are the 

eigenvalues of the matrix p{ay ® ay)p*{<Ty ® <7y), with 
Ai > Xi for i = 2,3,4, while -Efcp] (i, ^'o) is obtained 
as a sum of positive numbers, the latter reproduces the 
modulus of K{t) = A(p(t)). Consequently, even c^,[t) as 
derived from an individual quantum trajectory vanishes 
when A(t) changes sign at t = t^; this is associated with a 
jump of the phase 9{t) — arg(c^-!/)^]^) by tt. Thus, the dis- 
entanglement time ts is unambiguously encoded in each 
individual quantum state trajectory, as shown in Fig. [T] 

Let us finally describe an experimentally realizable 
setup for the direct monitoring of the open system en- 
tanglement evolution. For the optimal unraveling ([7|, 
the measurement record ([5| has the explicit form 



(2) 



Y.dt 



/7(- 



.(1) 



.dt 



(11) 



where, for notational simplicity, we let the operators' 
superscripts specify the subspace on which they act. 
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This experimental proposal, together with the time 
evolution equation ([9| and the finite time disentan- 
glement detection with a single trajectory, convey the 
strength of diffusive quantum trajectories as a complete 
and efficient method for the characterization of entan- 
glement under incoherent dynamics. Extensions of these 
ideas to other type of environments, higher dimensional 
systems, and considerations of finite detection efficien- 
cies for experimental realizations are relevant issues to 
be addressed in future work. 
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FIG. 2: Experimental setup for direct monitoring of two-qubit 
entanglement evolution under incoherent dynamics. One of 
the photons emitted by the two qubits passes through a phase 
shifter (PS) with 9 = — 6'opt + 7r/2 before combining with 
the other photon in a 50-50 beam splitter (BS). Each of the 
outgoing modes is then subject to homodyne detection with 
a local oscillator. 



Such measurement can be realized in optical experi- 
ments: Photons emitted by qubits 1 and 2, corresponding 
to Y^f-'* y^a'"^\ respectively, enter the detection 
scheme shown in Fig. [2j The photon coming from the 
second qubit passes through a phase shifter (PS) with 
9 = — ^opt + 7r/2 before impinging on a 50-50 beam split- 
ter (BS) , where it is combined with the photon from qubit 
one. Each of the outgoing modes is then subject to homo- 
dyne detection with a local oscillator with its phase set 
equal to ^loc = —tt/2. At the detectors, this procedure 
yields the real homodyne photocurrents [21 

h dt = \h-{-e''''^'a^+^ + a^'^ + h.c.),di + dCi , 

(12) 

h dt = ~iJ^-{e'''°-'af + a''^^ - h.c.)^ dt + dC2 ■ 

Here, dd^i are real, independent increments, correspond- 
ing to the detectors' shot noise. The complex cur- 
rents ( |Tl| ) are recovered from these real currents through 
Yi dt = (Ji dt - ih dt)/^/2 and Y2 dt = -e*^°p'(/i dt + 
i/2 dt)/^/2 . The conditional evolution of the state is re- 
constructed by substituting these, with the help of 
into ([3|. This then allows for an immediate evaluation 
of the pure state concurrence. The mixed state entangle- 
ment evolution of the system is recovered after averaging 
over the ensemble of entanglement records generated in 
this way. 
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